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Abstract. - Magnetic field induced antiferromagnetic phase of the underdoped cuprates is studied 
within the t — t' — J model. A magnetic field suppresses the pairing amplitude, which in turn 
may induce antiferromagnetism. We apply our theory to interpret the recently reported quantum 
oscillations in high magnetic field in ortho-II YBa2Cu306.5 [1] and propose that the total hole 
density abstracted from the oscillation period is reduced by 50% due to doubling of the unit cell. 
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Recently, Doiron-Leyraud et al. reported quantum 
oscillations in the clean underdoped cuprate ortho-II 
YBa2Cu306.5 in a high magnetic field at low tempera- 
tures [1] which imply the presence of small Fermi surface 
(FS) pockets. Their finding has shed new light on the 
nature of the anomalous spin gap phase of underdoped 
cuprates, a key to understanding the mechanism of high 
Tc superconductivity. The existence and topology of a 
FS [2,3] is one of the most important issues for underdoped 
cuprates. Assuming the normal state is paramagnetic as 
in zero field, angle resolved photoemission spectra on other 
underdoped cuprates [4] imply a FS consisting of four hole 
pockets centered on the Brillouin Zone (BZ) diagonals. 

In this Letter we propose that the reported oscillations 
' occur in an antiferromagnetic (AF) ordered state, hence 
the oscillation period corresponds to a Fermi pocket in 
the extended BZ. Our proposal is substantiated by ex- 
amining the interplay between AF and SCinat — i' — J 
model using a renormalized mean field theory (RMFT) [5] . 
A high magnetic field suppresses the pairing amplitude, 
which may in turn induce AF ordering. The FS prop- 
erty is studied using a phenomenological Green's function 
approach. The measured oscillation period in the pres- 
ence of AF ordering implies a total hole concentration of 
0.075 instead of 0.15 as interpreted by Doiron-Leyraud et 
al. [1]. The latter substantially exceeds the value of 0.10 
obtained from an empirical relation for Tc — of this 



compound [6] and a theoretical estimate of 0.04 based on 
band structure calculations [7]. The well ordered nature 
of this compound should enhance T^, and the empirical 
formula should overestimate, not underestimate, the hole 
concentration. We believe that our proposal gives a more 
reasonable hole density for this compound. 

The interplay of AF and SC order in underdoped 
cuprates has been a central issue [8-12]. Recent nuclear 
magnetic resonance data on multilayer cuprates have es- 
tablished that uniform AF order can coexist with SC order 
in the same Cu02 plane even in the absence of a magnetic 
field and have led Mukuda et al. to propose a revised phase 
diagram of the cuprates with coexisting AF and SC order 
in the underdoped region [13]. Generally in cuprates AF 
and SC orderings are separated by a spin glass region pre- 
sumably due to their crystalline disorder. Coexistence of 
uniform AF and SC order has been predicted from Varia- 
tional Monte Carlo (VMC) calculations on Gutzwiller pro- 
jected wavefunctions [10, 14]. We note that in the case of 
La2-xSrxCu04 neutron scattering data in an applied mag- 
netic field revealed subgaps in spin excitations [15], which 
were interpreted as evidence for a field induced AF or- 
dering in the SC state [16]. Our microscopic theory is 
consistent with their phenomenological analyses [16]. 

We consider the t — t' ~ J model on a square lattice 



p-1 



Wei-Qiang Chenl, Kai-Yu Yangl, T. M. Ricel,2, F. C. Zhangl 



[17,18], 

H = -Y, tij {clcj, + h.c.) + J ^ Si • S,- . (1) 

The hopping integrals are tij = t for the nearest neigh- 
bor (n.n.) pairs and tij = t' = —t/3 for the next n.n. 
pairs. Si is a spin-1/2 operator, and the spin exchange 
term (</ = t/3) connects n.n. pairs. The constraint of no 
double occupation at any lattice site is implied. We use 
a variational wavefunction [9,10] to study the interplay of 
AF (or spin density wave (SDW)) and SC ordering, 

I*) = Pzjjfo), |*o) = n ("k. + t^k«4.T^-k.i)|0), 

k,s=± 

(2) 

where the product of k runs over the reduced BZ, and 

Pjj = rii(-'^~^iT'^«i) Gutzwillcr projection operator. 

dksa is an annihilation operator of the SDW quasiparticles, 
with s = ± for the upper or lower SDW bands, 

rfk+.cr = COS (0k/2)ckCT - crsin (0k/2)ck+Qcr, (3) 
dk-,a = sin((/)u/2)ckff + (TCOs((/)k/2)ck+Q<T, 

where Q = (tt, tt) is the AF wave vector. 

To carry out the variation procedure we apply the 
Gutzwiller approximation to replace the effect of the pro- 
jection operator by a set of renormalization factors, which 
are determined by statistical counting [5, 19] and later 
improvement [11,20]. This approximation or the renor- 
malized mean field theory [5] incorporates the resonating 
valence bond physics proposed by Anderson [17] and com- 
pares well with the VMC results [19]. Furthermore, as has 
been emphasized recently by Anderson et ai, the RMFT 
describes many key features of the phase diagram [19] . Let 
(O) be the expectation value of operator O in the pro- 
jected state and (O)o be that in the corresponding 
unprojected state |^'o), then 

(clcja) = 5tS(cLc,.)o, {SJSJ) = 9l{S:SJ)o, (4) 

where t = x,y, z, and g's are the renormalization factors. 
In the AF state, electron populations in the two sublattices 
A (with net spin up) and B (with net spin down) are 
different. For the n.n. pair (ij), = 2(5(1 — S)/{1 — 
+ 4mo). For the next n.n. pairs on sublattice A or B, 
= g^_„ = gt{l + 6 + 2amo)/{l + 5- 2(jmo). In the 

above expressions, 5 is the hole density, mo = (sf (—1)'^ 
is the staggered magnetic moment in unit oi Hb in the 
state To address the interplay between SC and AF 

ordering, we use an improved Gutzwiller approximation of 
Ogata and Himeda for the spin exchange renormalization 
factors, which takes the effect of intersite correlations into 
account [11]. In their formalism, are also functions 
of mo and the pairing amplitude A in Eq. (6), and are 
given by 5^^^^ = gj/ 6'^{l+r]i)-'' , and = 5^(1-^772) with 



rj2 both small positive numbers, given in Ref. [21]. The 
difference between gl and g^'"^^ is due to the asymmetry 
in spin space. 

Within the Gutzwiller approximation, varying a pro- 
jected state for H in Eq. (1) reduces to varying an unpro- 
jected state l^'o) for a renormalized Hq, 

^0 = - E a^^^o (4cia + h.c.) + J ^ glSJS]. (5) 

We introduce two mean fields, Xx = Xy = X for hopping 
and i^x = ~^y = •A for d-wave pairing, 

X = E (^-FxaCicr)^ , Aa, = {Ci+x^Cil - Ci+xiCi|)o (6) 

The singlet SC order parameter lS.sc ~ .9t'A, with gj' = 
(5^1 + 5^1- )/2- The pairing amplitude and the SDW state 
described below define We choose a standard SDW 
form, 

cos (.^k) = ek/Ck, Ck = V4 + ^m,k' (7) 

where ek = —{2tgt + 3J^sX/4)7-i-(k) is the kinetic energy 
contribution of n.n. hopping including a self-energy term 
of X, with 5, = {gl + 25f 2')/3, and A„,k = m„ 2t'{g{^^ - 
g^lWk- fny is a variational parameter which determines 

™o, 7±(k) = cosfcj; ± cosfcj,, and 6'k = cos fcj, cos fc^ . The 
d-wave pairing variational functions are chosen as 

vl^ = (1 - Ck±/£;k±)/2, (8) 
Mk±Wk± = A^,(k)/2£;k±, 
a± = ±Ck - 4t'5t.^k - 
A, (k) = (3/4)J5,A„o7-(k), 

with £'k± = y^^ki + At.'(k)^, and A^jo is a variational pa- 
rameter. Note that because of the mean field dependence 
of the ^(-factors here, A^o 7^ A in general. 

The variational energy per site E = (H) /Ns can then 
be expressed in terms of the mean field parameters, 

E = -4gttx-St'grXi+At'{g^^^g^^)x2 (9) 
-{3gs/4)J{A''+x^)-2g,Jml 

where xi =^ -(1/2A?-,) Ek,s ^ks^k/^^ka, X2 = 

-(l/2A^s)Eks-sA™,k^ks6'k/CkSk,s- The variational 
ground state is obtained by minimizing E with respect 
to rriy, AyQ, X, A*, under constraint of the constant total 
hole concentration S. In Fig. 1, we plot the pairing mean 
field A, the SC order parameter Asc, and the staggered 
magnetization m = ^/gfmo in the projected state ]^), as 
functions of doping 5. At S < Sc = O-l, AF order co-exists 
with SC order; at 6 > 6c, the ground state has only SC 
order with m = 0. This result is similar to the variational 
Monte Carlo calculation which gives Sc — 0.10 for the 
t-J model [14]. 
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We now consider the effect of an external magnetic field 
on the phase diagram. Cuprates arc type II supercon- 
ductors, and a magnetic field penetrates in the form of 
vortices. The SC order outside the vortex cores is sub- 
stantially suppressed if the applied field Happ is compara- 
ble with the critical field 11^2 as in the experiments in Ref. 
[1] , but quantum oscillations may still occur [22] 

We use the RMFT to show that the suppression of the 
SC pairing amplitude due to the magnetic field may induce 
AF ordering. Consider the case 5 > 6c, so that the ground 
state at H, 
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is a pure SC state with the values in Eq. 
(6) Ao and mo = 0. We approximate the effect of Happ on 
the background SC region by a suppression of the pairing 
amplitude from Aq to an average value A = aAo, with A 
estimated below. Then Eq. (9) allows us to carry out the 
variation for a given A or equivalently for a fixed value 
of the variational parameter Avq. We have found that at 

5 > Sc. & suppression of Asc and hence A will induce AF 
ordering m 7^ 0. In Fig. 2 we plot m as a function of a for 

6 = 1.26c — 0.12. A may be estimated as follows. Let S^coh 
be the SC coherence length, and d the average distance of 
two neighboring vortices. We have d/S^coh ~ \/ Happ/Hc2, 
and A = 7r~^(d/2)-2 /^^^^^ ^ (r) rfr, where A (r) is the 
shape of the pairing mean field around a vortex. For 
a rough estimate, we use for s-wave the form, A (r) = 
Aotanh(r/^) [23]. The value of Hc2 of underdoped high 
Tc superconductors is about 160r [24], so the suppression 
ratio a at Happ = lOT (~ i?c2/16), 1ST (~ ifc2/9), and 
40T (~ i?c2/4) are 0.82, 0.72, and 0.56, respectively. The 
points corresponding to these values are denoted in Fig. 2 
with open circles. 

We turn now to the implications of our calculations for 
the interpretation of the recent experiments by Doiron- 
Leyraud et al. [1]. In their experiments, the Hall resistance 
of underdoped YBa2Cu306.5, as a function of magnetic 
field, shows clear oscillation in the vortex liquid phase. 
The FS area thus obtained is found by using Onsager re- 
lation and corresponds to small FS pockets with a charge 
carrier density of 0.0375. The compound is SC without AF 
ordering at low temperature in the absence of magnetic 
field. Estimates of the hole doping from the SC transi- 
tion temperature give a value 6 = 0.1. Since Fermi arcs 
have been observed in angle resolved photoemission spec- 
tra around the nodal points in Na-CCOC at similar hole 
doping [4], Doiron-Leyraud et al. interpreted the full FS 
to be composed of 4 such Fermi pockets. Here we point 
out that while the low temperature phase at Happ = 
is SC, the high magnetic field may induce AF ordering 
as discussed above. If this is the case, the area enclosed 
by the full FS will be double the pocket area. Thus the 
quantum oscillation may actually indicate a hole doping 
of only 0.075 instead of 0.15. 

Since quantum oscillations in the vortex liquid phase 
are known to reproduce those of the underlying normal 

state [22], we examine the FS pockets in the non-SC state. 
To this end we extend the phenomenological theory we 



developed earlier for the pseudogap state [25] to include 
AF ordering. Wc make an ansatz for the coherent part of 
the Green function (GF) which represents a quasiparticle. 



L0~ ^ks- Ss(k,cj) ' 



(10) 



where z is a numerical factor for the weight of the coherent 
part. The self energy takes the form 



E,(k,u;)=AV(k)/(u; + sCk), 



(11) 



where Cks and ^ks are given in Eqs. (7) and (8) with the 
replacement of m„ by toq. The Luttinger sum rule [26,27] 
for the number of electrons can be formulated as 



Afe = ^e[i?eG^^(k,a; = 0)], 



(12) 



ks 



where Q{x) is the Heaviside step function. Note that the 
boundary of the reduced BZ in the SDW state coincides 
with the Luttinger surface of zeroes of the corresponding 
GF in the absence of AF ordering. The FS of G\p depends 
on the staggered magnetization rriQ. In Fig. 3, we show 
the FS for several values of mo at a hole doping 6 = 0.075. 
In that case, the upper SDW band is completely empty 
and the lower SDW band is partially filled. A is chosen 
to be the value defined as in Eq.(6) that gives the lowest 
energy in the RMFT for a given mo- In such case the 
orbits of a quasiparticle in a magnetic field are confined to 
the reduced BZ and consist of two ellipses each composed 
of two half ellipses centered on opposite diagonals BZ. The 
period of quantum oscillations corresponds to one half of 
the total hole density. In the absence of AF order, i.e. 
mo = 0, the phenomenological GF has now four separate 
pockets near the nodal points and associating the period 
of the quantum oscillations with the area of a pocket leads 
to a total hole density twice as large. 

The electron GF can be obtained from the SDW GF. 
The resulting electron spectra weights relevant for angle 
resolved photoemission spectra measurements are shown 
in Fig. 3 by the width on the Fermi pockets. The spectral 
weight distribution inside the reduced zone at small mo 
is similar to our phenomenological theory for the pseudo- 
gap state [25] , i.e. a substantial weight only on the inner 
edge. The spectral weight outside the reduced BZ is very 
small but grows as mo increases. With increasing mo, the 
pockets become more ellipsoidal. 

In summary, we have used the renormalized mean field 
theory to show that a high magnetic field may induce 
antiferromagnetism in underdoped superconducting state. 
Applying our theory to the recently observed quantum os- 
cillations in a magnetic field in underdoped cuprates [1], 
we propose that the period measured in oscillations im- 
plies a hole density of 0.075, instead of 0.15 as originally 
proposed. 

Since the original manuscript was prepared, quan- 
tum oscillations have been observed in Y Ba2CuiOs [28]. 
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Again assuming 2 rather than 4 pockets gives a more rea- 
sonable vahio of the hole density 6 = 0.1. Lee has drawn 
our attention to his proposal of a doubled unit cell due to a 
staggered flux phase [29] as an alternative to the AF phase 
considered here. We have carried out a similar analysis for 
this case and found that it has a slightly lower energy than 
the AF phase. It may well be possible to observe the bro- 
ken symmetry in moderate fields (eg see Fig. 2) and so 
check for both possibilities experimentally. Renormalized 
mean field theory show that a high magnetic field may 
induce antiferromagnetism in undcrdopcd superconduct- 
ing state. Applying our theory to the recently observed 
quantum oscillations in a magnetic field in underdoped 
cuprates [1], we propose that the period measured in os- 
cillations implies a hole density of 0.075, instead of 0.15 
as originally proposed. 

After we submitted the paper, LeBoeuf et al. [30] pub- 
lished their new data which strongly suggest that at high 
magnetic fields in both YBCO compounds an unpaired 
normal state is stabilized with a longer period AF or other 
superlattice rather than the simple AF or flux ordering 
that we consider in this paper. This implies that the over- 
all negative sign of the Hall signal cannot be attributed 
to a contribution from vortex motion and is due rather 
to electron than to hole pockets. Millis and Norman have 
found that such an electron pocket may appear and domi- 
nate magnetotransport in the presence of an AF superlat- 
tice with an 8-fold period perpendicular to the chains at 
1/8 doping. [31] Further calculations extending our results 
are underway and will be reported elsewhere. 

* * * 

We thank L. Taillefer for providing us with a copy of 
their paper prior to publication and we thank him and P. 
A. Lee for discussions. This work was in part supported by 
RGC grant of HKSAR. FCZ wishes to thank Joint Theory 
Institute of University of Chicago and Argonne National 
Laboratory for its hospitality, where part of the work was 
carried out. Support by the MANEP program of the Swiss 
National Foundation and a Visiting Professorship at the 
University of Hong Kong is also gratefully acknowledged. 

REFERENCES 

[1] Doiron-Leyraud N. et al., Nature, 447 (2007) 565. 
[2] Ding H. et al.. Nature, 382 (1996) 51. 
[3] Marshall D. S. etal, Phys. Rev. Lett, 76 (1996) 4841. 
[4] SnEN K. M. et al, Science, 307 (2005) 901. 
[5] Zhang F. C. et al, Supercond. Sci. Technol, 1 (1988) 
36. 

[6] Liang R. X. et al, Phys. Rev. B, 73 (2006) 180505(R). 
[7] Bascones E. et al, Phys. Rev. B, 71 (2005) 012505. 
[8] For recent reviews see Birgeneau R. J. et al, J. 

Phys. Soc. Jpn., 75 (2006) 111003; Tranquada J. M., 

arXiv, cond-mat/0512.115 (2005) . 
[9] Lee T. K. and Feng Shiping, Phys. Rev. B, 38 (1988) 

11809. 



[10] Chen G. J. et al, Phys. Rev. B, 42 (1990) 2662. 
[11] Ogata M. and Himeda A., J. Phys. Soc. Jpn., 72 (2003) 
374. 

[12] Zhang S. C, Science, 275 (1997) 1089. 
[13] MUKUDA H. et al, Phys. Rev. Lett., 96 (2006) 087001. 
[14] Himeda A. and Ogata M., Phys. Rev. B, 60 (1999) 
R9935. 

[15] Lake B. et al. Science, 291 (2001) 1759. 
[16] Demler E. et al, Phys. Rev. Lett, 87 (2001) 067202. 
[17] Anderson P. W., Science, 235 (1987) 1196. 
[18] Zhang F. C. and Rice T. M., Phys. Rev. B, 37 (1988) 
3759. 

[19] Anderson P. W. et al, J Phys. Condens. Matter, 16 
(2004) R755. 

[20] Can J. Y. et al, Phys. Rev. Lett, 94 (2005) 067005. 

[21] The small parameters 771 and 772 in (/s and are given 
by??i =2X/(l-52+4mg)2,?72 = 8ml{2Xi+Xf)/{8ml + 
X2), X = 6^ (A^ - x') + 4m§X2 + Xl/4, X2 = + A^ 
and Xi = 3X2(1 - <5)V[(1 - <5' + 4mg)(l + ??i)3]. 

[22] Maniv T. et al., Rev. Mod. Phys., 73 (2001) 867. 

[23] Gygi F. and SCHLijTER M., Phys. Rev. B, 43 (1991) 
7609. 

[24] Ong N. p. and Wang Yayu, Physica C, 408 (2004) 11. 

[25] Yang K. Y. et al, Phys. Rev. B, 73 (2006) 174501. 

[26] Konik R. M. et al, Phys. Rev. Lett., 96 (2006) 086407. 

[27] Abrikosov a. a., Gorkov L. P., and Dzyaloshinaskii 
I. E., Methods of Quantum Field Theory in Statistical 
Physics, revised ed., edited by Silverman R. A. (Dover, 
New York) 1975; Dzyaloshinskii I., Phys. Rev. B, 68 
(2003) 085113. 

[28] Yelland E. a. et al, Phys. Rev. Lett, 100 (2008) 
047003; Bangura A. F. et al, Phys. Rev. Lett, 100 
(2008) 047004. 

[29] Lee P. A., arXiv, cond-mat/0708.2115 (2007) . 

[30] LeBoeuf D. et al. Nature, 450 (2007) 533. 

[31] A. J. Millis and M. Norman, Phys. Rev. B, 76 (2007) 
220503(R). 



0.35 



0.3 -\ 




5 

Fig. 1: Pairing mean field A (dashed line), superconducting 
order parameter Age and antiferromagnetic staggered magne- 
tization m (solid lines) as functions of hole concentration 5 for 
t-t' - J model, with J = t/3, t' /t = -1/3. 
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Fig. 2: Staggered magnetization m obtained from Eq. (9) of 
RMFT as a function of the ratio of the suppressed pairing 
amplitude and the zero-field pairing amplitude a = A/Ao at 
doping 6 = 1.2Sc_{Sc = 0.10) (t = 0.3 eV, J = 0.1 eV, t' = 
—0.1 eV), where A is the spatially averaged pairing mean field 
A (r) as discussed in context, Ao is the pairing amplitude at 
zero field. Circles indicate the estimated corresponding applied 
fields (see the text). 
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Fig. 3: Evolution of Fermi surface in one quadrant of BZ cal- 
culated from Eq. (10-12) at hole doping S — 0.075. The period 
of quantum oscillations is given by the area of the closed Fermi 
pocket, which is (5/4 for (a) in the absence of AF order and is 
5/2 for (b)-(d) in the AF ordered state. The thickness of the 
curves represents the spectral weight. 
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